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Abstract. We present some recent existence results for the varia- 
tional model of crack growth in brittle materials proposed by Franc- 
fort and Marigo in 1998. These results, obtained in collaboration 
with Francfort and Toader, cover the case of arbitrary space dimen- 
sion with a general quasiconvex bulk energy and with prescribed 
boundary deformations and applied loads. 
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Introduction 

In 1998 Francfort and Marigo [12] introduced a variational model for 
the quasistatic growth of brittle cracks in elastic materials. This model is 
based on Griffith's idea [13] that at each time the equilibrium of a crack 
is determined by the balance between the elastic energy released when 
the crack grows and the energy dissipated to produce a new portion of 
crack. This model determines the crack path on the basis of an energy 
criterion, and describes also the process of crack initiation. 

1. The model 

The reference configuration. The reference configuration is a bound- 
ed open set of W 1 with Lipschitz boundary dfL We fix a partition 
d£l = dci£l U djsffl. On the Dirichlet part dc>£l of the boundary we 
prescribe a time-dependent boundary deformation, while the Neumann 
part is traction free. 
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The crack. In this model a crack is any a countably (TC a ~ l ,n — 1)- 
rectifiable set T contained in 17 and with H n_1 (r) < +oo, where 'Yi n ~ 1 
is the (n — 1) -dimensional Hausdorff measure (see [10] for the definitions 
of these notions of geometric measure theory). No a priori assumption 
is made on the shape or on the topology of the crack. 

The energy dissipated to produce the crack V depends on the inter- 
atomic bonds broken in this process. We assume that it can be written 
as 



where v is a unit normal vector field on T. The function k{x,v) is 
called the toughness of the material. We assume that it is even in v 
and satisfies the standard hypotheses which guarantee the lower semi- 
continuity of /C. Since k(x,u) depends both on the position x and 
on the orientation z/, this model covers heterogeneous and anisotropic 
materials. 

The deformation. Under the effect of the body forces and of the 
boundary conditions the uncracked part Q \ T of the body undergoes a 
deformation, described by a function u: Q, \ T — > R n . 

To give a precise mathematical formulation of the problem, we have 
to choose a suitable function space for the deformations u. Since T 
is not prescribed, but has to be determined on the basis of an energy 
criterion, it is not convenient to work in spaces of functions defined on 
domains depending on T . 

It is more convenient to consider u as a function which is defined 
almost everywhere on Q and is discontinuous on an (n — 1) -dimensional 
set. Spaces of discontinuous functions of this kind have been introduced 
by De Giorgi and Ambrosio (see [9] , [1] , and [4] ) to study a large class of 
free discontinuity problems, where one minimizes functionals depending 
on the discontinuity set of the unknown function. 

For the deformations we use here the space GSBV P {Q; W 1 ) defined 
as the space of all functions u G GSBV(Q;R n ) with H n_1 (J(u)) < +oo 
and Vu G Z7(ft;]R nxn ) (see [1] or [4] for the definitions). 

Without entering into details, for the purposes of this exposition it is 
enough to know that for every u G GSBV P {£1; W 1 ) we can define, in a 
mathematically precise way, 

• the jump set J{u) , which is a countably {ri n ~ l , n — l)-rectifiable 
subset of n with U n ~ l {J{u)) < +oo; 
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• the approximate gradient Vu, which belongs to L p (U;W nxn ); it 
coincides a.e. in £1 with the ordinary gradient of u, when J(u) 
is closed in and u is smooth in fi \ J(u) ; 

• the trace of u on dil , defined through the notion of approximate 
limit. 

The bulk energy. We assume that the material is hyperelastic in the 
uncracked part Q\T of the body. For a deformation u G GSBV P (£1; R n ) 
with J(u) C T , the elastic energy stored in Q \ T is given by 

W(Vu) := / VF(x,Vn(x))dx, 
jn\r 

where W(x,£) is a given function depending on the material. To guar- 
antee lower semicontinuity, we make the usual assumption that £ i— > 
W(x,£) is quasiconvex and satisfies the standard coerciveness and growth 
conditions of order p > 1 . Unfortunately these conditions do not include 
the case of finite elasticity, where W(x,£) = +oo when det(£) < 0. The 
extension of our results to this case remains an open problem. 
It is useful to consider also the functional 

: = f W(x,$(x))dx, (1.1) 
defined for every vector field $ G L p (Q;R nxn ) . 

The body forces. For every time t > the applied load is given by 
a system of t-dependent body forces. We assume that these forces are 
conservative and that their work on the deformation u is given by 

F(t)(u) := / F(t,x,u(x))dx, (1.2) 
Jn\r 

where F satisfies suitable regularity and growth conditions. In partic- 
ular we assume that there exist an exponent q > 1 and two constants 
a > and (3 > such that 

-F(t,x,u) > a\u\ q -j3. (1.3) 

This condition, which ensures the coeciveness of the total energy that 
we will introduce in (1.4), says, in a very weak sense, that the forces 
are attractive at large distances. The reason for this assumption is 
that we want to prove an existence result for the quasistatic evolution 
for arbitrarily large times. During the process of crack growth it may 
happen that the crack breaks the body into serveral pieces, some of 
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which may be disconnected from the Dirichlet part of the boundary. 
These pieces run the risk of being sent to infinity by the body forces, 
unless a condition like (1.3) is satisfied. 

The boundary deformations. For every time t > we prescribe a 
boundary deformation w{t) on \ T(i) , where T(t) is the unknown 
crack at time t. Since the deformation may be discontinuous on the 
crack, it makes no sense to prescribe the boundary deformation also on 

d D nnr(t). 

We assume that w{t) is sufficiently regular on dr>£l, so that the pres- 
ence of the crack is not imposed by prescribing a "strong discontinuity" 
at the boundary. Moreover, we also assume that t i— > w(t) is sufficiently 
regular with respect to time. 

The admissible configurations and their stability. Given a bound- 
ary datum w: dpti — > K, the set A(w) of all admissible configurations 
with boundary deformation w is defined as the set of all pairs defor- 
mation-crack (u, T), where the crack T is a countably (7Y n_1 ,n — 1)- 
rectifiable set contained in Q, with W ra_1 (r) < +oo, while the deforma- 
tion u belongs to GSBV P (U; W 1 ) nL'fOjK") and satisfies J(u) C T 
and u = w on do^ \ T . 

The total energy at time t of a configuration (u,T) £ A(w{t)) is given 

by 

£(t){u, T) := W(Vu) + /C(r) - ^(t)(u) , (1.4) 

the sum of the stored elastic energy of the deformation u , of the energy 
dissipated to produce the crack T , and of the opposite of the work done 
by the body forces. 

In the spirit of Griffith's original theory, a configuration (u(t),T(t)) 
is globally stable at time t if (u(t),T(t)) G A(w(t)) and 

£(t)(u(t),T(t))<£(t)(u,T) (1.5) 

for every (u,T) <G A(w(t)) with T D T(t). In other words, the energy of 
(■u(i),r(t)) can not be reduced by choosing a larger crack and, possibly, 
a new deformation with the same boundary condition w(t) . 

We observe that, from the point of view of mechanics, it would be 
preferable to consider a local version of the minimality condition (1.5). 
Indeed, a configuration is still in equilibrium if inequality (1.5) is sat- 
isfied only when T is close to F(t) and u is close to u(t). So far the 
mathematical theory has been fully developed only for globally stable 
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configurations, and only partial results have been obtained for other 
local equilibria (see [8]). 

Irreversible quasistatic evolution. An irreversible quasistatic evo- 
lution of globally stable configurations is a function t t— > (u(t),T(t)) 
defined for t > which satisfies the following conditions: 

(a) global stability: for every t > the pair (u(t),T(t)) is globally 
stable at time t; i.e., (u(t),T(t)) G A(w(t)) and 

£(t)(u(t),T(t)) <£(t)(u,r) 

for every (u,T) G A(w(t)) with T D r(i); 

(b) irreversibility: T(ti) C T(t 2 ) for < ii < i 2 ; 

(c) energy balance: in every time interval the increment of the stored 
elastic energy plus the energy dissipated to extend the crack is 
equal to the work done by the external forces, including the 
unknown forces acting on dpti \ r(i) . 

In (c) the increment of the stored energy in the interval [ii, <2] is given 
by W(Vu(t2)) — W(Vu(ti)), while the energy dissipated by the crack 
growth is /C(r(t2) \ r(ii)). Using the Euler equation for the global 
stability condition it is possible to show that the work of all external 
forces involved in the problem is given by 



where dots denote time derivatives, while dW and dJ-(s) are the differ- 
entials of the functionals defined by (1.1) and (1.2) on L p (Q;R nxn ) and 
L q (Q; R n ) , respectively. It turns out that (c) can be written equivalently 
in the following differential form: 

(c) energy balance: the function t i— > E(t) := £(t)(u(t),T(t)) is lo- 
cally absolutely continuous and its time derivative E(t) satisfies 



E(t) = (dW(Vu(t)),Vw(t)) - {<LF(t)(u(t)),w(t)) -.F(t)(u(t)) 



for a.e. t > 0. 

There is strong numerical evidence that, while the function t >— > 
£(t)(u(t),F(t)) is continuous, the functions t t— > W(Vu(£)) and t \— > 
/C(r(i)) may present some jump discontinuity. The smallness of the set 




{(dW(V«(s)), Vw(a)> - (dF(a)(«(a)),w(a))}da + 
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of discontinuity times is still an open problem, as well as the regularity- 
properties of these functions out of the discontinuity set. 

2. The existence results 

The results. The most general existence result is given by the following 
theorem, proved in collaboration with Francfort and Toader [6]. 

Theorem 2.1. Let (no, To) be a globally stable configuration at time 
t = . Then there exists an irreversible quasistatic evolution of globally 
stable configurations t i— > (u(t),T(t)) with (u(0),r(0)) = (uo,Tq) . 

Previous results on this subject have been obtained in [7] in the case 
n = 2 for a scalar-valued u and for W(£) = |£| 2 , which corresponds 
to the antiplane case in linear elasticity. In that paper the admissible 
cracks are assumed to be connected, or to have a uniform bound on 
the number of connected components. This restriction, which has no 
mechanical justification, allows to simplify the mathematical formulation 
of the problem. These results were extended to the case of planar linear 
elasticity by Chambolle [5]. 

A remarkable improvement, still in the case of a scalar- valued u, was 
obtained by Francfort and Larsen [11], who were able to remove the 
restrictions on the number of connected components of T and on the 
dimension of the space. 

Theorem 2.1 is the first result which covers the case of a vector- valued 
u in any space dimension. Moreover, it includes time-dependent body 
forces, that could not be treated by the methods used in the previous 
papers. 

The incremental problems. As in the other works on this subject 
[7], [5], [11], the proof of Theorem 2.1 is obtained by time discretization. 
For every k we fix an increasing sequence (t l k )i>o , with 

o = 4 < t\ < ■ ■ ■ < f- 1 < 4 - +cx) . 

We also assume that 

lim sup(4-4- 1 ) = 0. 

fe^OO j 

The approximate solutions (u l k ,r l k ) at time t\ are defined by in- 
duction: we set (u°,r°) := (uq,Tq), and, assuming that {u 1 ^ 1 1 T l k ~ l ) 
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is given, we define (u k ,F k ) as a solution to the incremental minimum 
problem 

min {£(tl)(u, r) : (u, T) € A{w(t{)), T D T^ 1 } . 

Note that the external forces and the boundary condition refer to the 
updated time t\, while the constraint T D T^ -1 refers to the previous 
time t l k ~ . The existence of a solution can be easily deduced from the 
GSBV compactness theorem of [2]. 

Passing to the limit. To pass from the discrete-time formulation of 
the incremental problems to the continuous-time formulation of the 
quasistatic evolution, we consider the piecewise constant interpolations 
u k (t) and Y k (t) defined by 

u k {t):=u%, r fc (*):=r* fc for t* fc < * < 4 +1 . 

In the first step of the proof we exploit suitable compactness properties 
and pass to the limit as k — > oo along a subsequence: 

u k {t) u(t) , r k (t) -> r(t) . 

For the sequence of deformations u k {t) we use the following notion of 
convergence related to the space GSBV P (Q; M n ) : 

u k(t) ~^ u {t) pointwise a.e. on £1 , 
Vu k {t) -± Vu{t) weakly in R nxn ) , 

limsup fc TL n ~ l {J (u k {t)) < +oo . 

The results of [2] and [3] provide good compactness and semicontinuity 
properties for this kind of convergence. 

For the sequence of cracks T k (t) the standard notion of convergence 
in the Hausdorff metric has good compactness properties, but it does 
not lead to the proof of the global stability of the limit pair (u(t),T(t)) , 
unless n = 2 and the sets Y k {t) are closed and have a uniformly bounded 
number of connected components. 

To treat the general case we have to introduce a new notion of conver- 
gence of sets, called a p -convergence, related to the notion of jump sets of 
SBV functions. For this notion of convergence we prove a compactness 
result similar to Helly's theorem for monotone functions, which allows 
to extract a subsequence, independent of t , still denoted T k , such that 
Tfc(t) a p -converges to T(t) for every t > 0. 
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Global stability and energy balance. It remains to prove that the 
limit pair (u(t),F(t)) satisfies the global stability condition (a), the ir- 
reversibility condition (b), and the energy balance (c). 

The global stability is proved by using a slight modification of a 
very deep approximation result proved by Francfort and Larsen [11] 
and known as Jump Transfer Theorem. 

Since, by construction, Tk(t\) C Tkfo) for t\ <t2, the same property 
holds for T(t) . Therefore the irreversibility condition is satisfied. 

In the proof of the energy balance there are two difficulties. The first 
one is the fact that, when we pass from the discrete-time formulation to 
the continuous-time formulation, we are led to approximate a Lebesgue 
integral by Riemann sums. 

The second difficulty is due to the fact that Vufc(i) converges to 
Vu(t) only weakly in L p (0,;M. nxn ) , and we have to pass to the limit 
in the nonlinear expression d^W(x,Vuk(t)) , where d(W denotes the 
nxn matrix of the partial derivatives of W(x,£) with respect to the 
components of £ 

The first problem is solved by using a suitable extension of a classical 
result which ensures that every Lebesgue integral can be approximated 
by carefully chosen Riemann sums. 

The second problem is solved by using an extension to GSBV p (£l; R n ) 
of the following lemma, which is new and interesting also in the Sobolev 
space W^tyR"). 

Lemma 2.2. Suppose that Uk — 1 u weakly in W 1,P (Q; R n ) and that 



Then d/:W{x,Vu k {x)) d^W(x,Vu(x)) weakly in i/ (Q; R nxn ) , with 
P' = P/{P ~ 1) • 

If W(x,£) is strictly convex in £, using [14, Theorem 2] one can prove 
that V«fc converges in measure to Vit , and the result follows easily from 
the growth conditions. We notice that in the quasiconvex case the result 
holds even if Vu k does not converge in measure to Vn. 
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